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ARO  Final  Report  -  Harvard  Lead  Phase 


Statement  of  the  problem  studied:  During  the  Harvard-lead  Phase,  the  ARO-IARPA  MQCO  team 
investigated  three  types  of  spin  qubits,  both  realized  in  coupled,  gate-defined  quantum  dots  in  high- 
mobility  GaAs  heterostructures.  The  first  is  the  Loss-DiVincenzo  (LD)  qubit,  where  a  single  electron  spin 
defines  the  logical  qubit.  The  second  is  the  singlet-triplet  (ST)  qubit,  where  the  two  zero-spin-projection 
subspace  of  two  spins  in  two  dots,  the  singlet  and  m=0  triplet — constitute  the  qubit.  The  third  is  the 
exchange-only  (EO)  qubit  comprising  three  spins,  where  exchange  between  electrons  1  and  2  gives  one 
rotation  axis,  and  exchange  between  2  and  3  gives  a  second  axis.  The  Bloch  spheres  for  the  three  qubits 
are  shown  in  Fig.l.  Measurements  are  carried  out  by  fast  charge  sensing  using  proximal  charge  detectors 
or  dispersive  charge  read-out.  Conversion  of  spin  states  to  charge  readout  takes  advantage  of  the  Pauli 
principle  preventing  double  occupancy  of  the  ground  orbital  state  of  a  dot.  Experiment  is  performed  in  a 
dilution  refrigerator  at  millikelvin  temperatures. 


Fig.  1.  Bloch  spheres  and  elementary  spin  representations  for  Loss-Divincenzo  (LD),  Singlet-triplet 
(ST)  and  Exchange  Only  (EO)  qubits  [from  left  to  right].  The  experimental  needs  for  qubit  control  are 
illustrated  below. 

The  reason  for  all  of  these  different  versions  is  that  the  development  of  the  qubit  itself  is  under 
development,  given  the  early  stage  of  the  technology.  ST  qubits  are  the  most  advanced,  EO  qubits  are  the 
least.  All  have  some  advantages  and  some  weaknesses. 

Pulses  for  initialization  and  manipulation  are  applied  to  electrostatic  gates  using  room-temperature 
arbitrary  waveform  generators  (AWG’s).  For  the  LD  qubit,  single  qubit  rotations  are  done  using  electric 
dipole  spin  resonance  (EDSR),  which  requires  pulsed  microwaves  at  frequencies  roughly  from  1  to  20 
GHz.  Charge  state  readout  uses  fast  amplifiers,  located  either  at  the  4K  stage  of  the  refrigerator  or  room 
temperature.  For  ST  and  EO  qubits,  ideally  square  voltage  pulses  on  gates  are  used.  Circuit  boards  for 
bring  GHz-scale  pulses  onto  the  chip  have  been  the  focus  of  considerable  effort  during  this  period.  An 
example  of  circuit  board  developed  during  the  Harvard-lead  period  is  shown  in  Fig.  3.  All  measurements 
are  carried  out  at  low  temperature  (mK).  A  new  cryogen-free  refrigerator  is  shown  in  Fig.  4. 
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Fig.  2.  Fast  readout  electronics  based  on  high-Q  resonant  circuit. 


Fig.  3.  High-bandwidth  circuit  boards  designed  for  spin  qubits. 
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Fig.  4.  A  cryogen-free  dilution  refrigerator  (opened). 
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Summary  of  Accomplishments 


Summary  of 
Accomplishments 


Measurement 


Results 


Notes 


Ti  Measurement 


Full  Control 


State  Tomography 


Universal  Entangling 
Gate  definition 


State  Tomography  and 
duration  of  gate 


See  Fig.  3.1.3 


See  Fig.  3.1.4 


See  Sec.  3. 


Tl  left  =  4.9  ±  1.7  ms 
Tl  right  =  3.8  ±  0.7  ms 


Poo 


Poi 


Pio 


Pll 


Detection  Time 


Detection  Fidelity  (ability 
to  distinguish  between  the 
two  qubit  states) 


See  Figs.  3.2.1-3.2.2 


In  progress 


P  left  0  =  0.950  ±  0.5  See  notes  in  Sec.  3.2 
P  left  1  =  0.780  ±  4 
P  right  0  =  0.957  ±  0.4 
P  right  1  =  0.777  ±  0.9 


0,95 


0,05 


0,78 


0,22 


~  2  ms  for  tunneling  event 


Upper  bound  for  the 
fidelity  after  a  quarter 
period  (sqrt  swap),  is  the 
fourth  root  of  0.977,  ~ 
0.994. 


Ti  Measurement 

1-10  ms  (See  Fig.  3.3.2) 

State  Tomography  and 
duration  for  gate 

Fidelity  0.97 

typical  gate  time  -  1 0  ns 

Process  Tomography  for 
(te/2)  gate 

Fidelity  95% 

3.4  Two  ST  Qubits 


See  Fig.  3.3.6 


See  Figs.  3.3.8  and  3.3.9 


See  Sec.  3.4 
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Summary  of 
Accomplishments 

Measurement 

Results 

Notes 

CPhase(7i) 

Include  Truth  Table 

ST  Readout 

See  Sec.  3.4 

Poo 

0,99 

See  Fig.  3.3.5 

Pol 

0,03 

See  Fig.  3.3.5 

PlO 

0,01 

See  Fig.  3.3.5 

Pll 

0,97 

See  Fig.  3.3.5 

Detection  Time 

0.8  jus 

Detection  Fidelity  (ability 
to  distinguish  between  the 
two  qubit  states) 

0.98 

3.5  Single  EO  Qubit 

See  Sec.  3.5 

Readout  Fidelty 

Singlet  Fidelity  0.96 

Triplet  Fidelity  0.90 

See  Fig.  3.5.2 

State  Tomography 

Average  fidelity  =  0.75 

See  Fig.  3.5.3 
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3.1.  Single  LD  Qubit 
T1  measurement 


Fig.  3.1.1  (A)  Charge  stability  diagram,  with  charge-sensing  current  Ilqpc  shown  in  color  scale  as  a 
function  of  voltages  applied  to  gates  LP  and  RP  (a  background  plane  has  been  subtracted).  The 
occupation  in  the  left  and  right  dots  is  indicated  by  numbers  in  brackets.  (Inset)  Scanning  electron 
micrograph  of  a  device  similar  to  the  one  used  in  our  experiment.  Gates  LP  and  RP  are  connected  to 
high-frequency  lines  via  bias-tees.  The  direction  of  Bext  is  indicated.  (B)  Electrochemical  potential 
diagrams  showing  the  double-dot  configuration  in  the  two  read-out  stages  [positions  (1)  and  (2)  in  (A)]. 
Tunnel  events  that  occur  for  a  H  state  are  indicated.  From  Ref.  [1]. 
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3.1.  Single  LD  Qubit 


Fig.  3.1.2  (Left)  Pulse  sequence  shown  on  charge  stability  diagram  including  the  (2,0)  charge  region. 
(Right)  Energy  diagram  showing  the  pulse  sequence.  From  Ref.  [1]. 


B  C 


Fig.  3.1.3  (A)  Single  dot  relaxation  as  a  function  of  wait  time.  Grey  lines  are  fits  to  p  ■  e~t/T1  +  a.  (B) 
Flistogram  showing  the  distribution  of  the  time  Td  it  takes  a  spin-down  electron  to  tunnel  out.  The  red 
line  is  an  exponential  fit  from  which  we  can  find  the  decay  rate.  Inset:  real  time  trace  of  one  of  the  read 
out  stages  indicating  the  detection  time  (Td)  and  the  event  time  (Te).  (C)  Flistogram  showing  the 
distribution  of  the  time  Te  it  takes  a  spin-up  electron  to  tunnel  back  into  the  empty  dot.  The  red  line  is 
an  exponential  fit  from  which  we  can  extract  the  decay  rate.  From  Supplemental  Material  to  Ref.  [1]. 
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3.1.  Single  LD  Qubit 
Full  Qubit  Control 


30  ns  90  ns 


Fig.  3.1.4  Ramsey  signal  as  a  function  of  the  relative  phase  between  the  two  rf  bursts  for  t  =  10ns 
(crosses)  and  150  ns  (circles).  Gray  dashed  line  is  a  best  fit  of  a  cosine  to  the  data.  From  [2]. 
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3.2  LD  Two  Qubit  Operation 

Entangling  operation:  Exchange  operation  on  two  LD  qubits. 
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Fig.  3.2.1  (A  to  D)  Two-qubit  exchange  gate  on  a  full  set  of  input  states.  The  four  panels  correspond  to 
four  different  mixtures  of  initial  states,  as  indicated,  taken  with  otherwise  identical  settings.  Again,  spin- 
down  injection  probabilities  are  below  50%.  Gray  lines  are  fits  to  damped  oscillations,  including  a 
correction  for  pulse  imperfections.  We  first  fit  PTj  in  (A)  and  P^t  in  (B)  and  allow  only  the  amplitude  and 
offset  of  the  oscillations  to  change  for  the  other  probabilities  in  the  respective  panel.  In  (C)  and  (D),  we 
use  the  fit  parameters  of  (A)  and  allow  only  amplitude  and  offset  to  change.  The  oscillations  in  (A)  and 
(B)  run  out  of  phase  with  each  other  for  longer  wait  times.  We  attribute  this  to  subtle  distortions  of  the 
pulses  arriving  at  the  sample  due  to  the  bias  tees  (22).  (E  and  F)  Visualized  theoretical  and 
experimental  truth  tables  for  a  tt  rotation  and  a  2tt  rotation  of  the  exchange  oscillation.  From  Ref.  [1]. 
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3.2  LD  Two  Qubit  Operation 

Entangling  operation:  Exchange  gate  with  single  qubit  operations 


(a) 


Avpl+n 


Initialize  (A) 

6=0 

Rotate  IB) 

W 

e»0 

SWAP  (C) 

e»0 

*2 

Rotate  IB)  Readout  (A) 


Single-qubit  gate 


Two-qubit  gate 


Fig.  3.2.2.  (a)  Scanning  electron  microscopy  image  of  the  device  fabricated  on  top  of  an  AIGaAs/GaAs 
heterostructure  showing  the  Ti=Au  gates  (light  gray)  and  the  split  cobalt  (Co)  magnet  (yellow) 
separated  from  the  gate  contacts  by  a  calixarene  layer.  Gates  R  (right)  and  L  (left)  control  Nr  and  Ni_;  C 
(center)  controls  the  interdot  tunnel  coupling  t.  Fast  voltage  pulses  are  applied  to  the  Co  and  PL  gates. 
A  MW  voltage  Vac  is  applied  to  the  upper  part  of  the  magnet.  Gqpc  is  measured  by  modulating  the  Pl 
gate  voltage  Vpl.  (b)  Cycle  of  the  two  qubit  gate  operation  with  source  (S),  drain  (D),  left  (L),  and  right 
(R)  QDs.  From  Ref.  [3].  (C)  Pulsing  voltages  to  demonstrate  the  quantum  operation  of  the  single-qubit 
and  two-  qubit  gate.  The  voltage  pulses  are  applied  to  switch  between:  A^B— >A  (controlled  single-spin 
rotations,  left  panel)  and  A— >B— >C— >B— >A  (quantum  circuit  comprising  controlled  single  spin  rotations 
and  exchange  operation,  right  panel).  The  control  cycle  for  the  controlled  single  spin  rotations  consists 
of  initialization  (A),  spin  rotation  (B)  and  readout  (A).  The  control  cycle  for  the  multiple-qubit  gate 
consists  of  initialization  (A),  spin  rotation  (0=3tt/2)  (B),  spin  exchange  or  "SWAP"  (C),  spin  rotation  (0 
=tt/2)  (B)  and  singlet  readout  (A). 
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3.3  Single  ST  Qubit 
T1  Relaxation 


Fig.  3.3.1  Capacitively  Coupled  two-ST  qubit  device.  From  Ref.  [4], 


Typical  Operating  Range 


Fig.  3.3.2  Qubit  relaxation  time  T1  as  a  function  of  detuning  into  the  measurement  position. 
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3.3  Single  ST  Qubit 
Definition  of  Fidelity 

In  order  to  quantitatively  interpret  sensor  values  for  state  tomography,  it  is  important  to  precisely 
determine  the  RF  sensor  response  that  corresponds  to  a  |  S)  or  a  |  T0)  state.  Because  the  state  preparation 
is  imperfect,  it  is  in  general  difficult  to  accurately  measure  these  values  experimentally.  To  provide  exact 
calibrations  for  |  S)  and  |  T0)  ?  we  exploit  the  fact  that  our  sensor  is  capable  of  single  shot  readout. 
Histograms  of  sensor  values  for  typical  measurements  yield  a  double-peaked  curve-  one  peak  corresponds 
to  \  T0)  and  one  to  |  S)  (Fig.  3.3.3).  In  order  to  calibrate  the  sensor  we  first  measure  1 1  at  the 
measurement  point  by  preparing  a  state  that  is  majority  |  T0)  (done  with  a  7i -pulse  around  the  x-axis)  and 
fitting  the  sensor  signal  to  a  decaying  exponential  function  of  time  elapsed  during  measurement  (Fig. 
3. 3. 3. a).  We  note  that  the  measured  value  of  Ti  is  a  strong  function  of  the  power  of  the  RF  excitation  used 
to  read  the  conductance  of  the  sensing  QD.  With  prior  knowledge  of  77,  we  use  a  procedure  similar  to  that 
described  in  Barthel  et.  al  (14)  to  optimize  the  measurement  time  given  our  signal  to  noise  ratios  and  77. 
This  process  is  repeated  several  times  per  day  to  check  for  drift.  We  recalibrate  the  sensor  signals  that 
correspond  to  |  S)  and  |  T0 )  for  each  dataset  (typically  10  minutes  of  acquisition  time).  For  each  set,  we 
prepare  a  histogram  of  all  observed  sensor  values.  The  presence  of  several  reference  measurements  in 
each  dataset  guarantees  that  there  will  be  a  significant  fraction  of  both  |  S )  and  |  T0)  .  We  then  fit  this 
double  peaked  curve  to  an  analytic  expression  corresponding  to  a  weighted  sum  of  two  Gaussians  with 
some  filling  in  due  to  77  decay  during  measurement  (Fig.  3.3.3.b,  purple  line)  as  in  ref.  14.  From  this,  we 
extract  the  expected  sensor  distribu|  T0)  tions  for  |  S)  and  |7"0)  (blue  and  red  lines  in  Fig.  Sib, 
respectively),  as  well  as  the  fractions  of  |  S )  and  |  T0 )  present.  The  centers  of  the  two  distributions 
correspond  to  the  sensor  signals  that  will  be  measured  for  pure  |  S )  and  pure  |  T0)  9  and  using  these  values 
we  can  accurately  scale  the  tomography  data.  We  note  that  this  procedure  is  insensitive  to  the  percentages 
of  |  S)  and  |7"0)  .  In  our  state  tomography  only  expectation  values  are  needed,  so  the  single-shot 
capability  of  our  readout  is  not  necessary  beyond  this  calibration.  Nonetheless,  we  note  that  for  the  data 
presented,  we  measured  readout  fidelities  of  97%  and  98%  for  the  left  and  right  qubits,  respectively. 


Fig.  3.3.3.  Singleshot  Readout:  a,  The  difference  in  sensor  signal  between  |  S )  and  |  T0)  is  fit  to  a 
decaying  exponential  to  determine  Ti,  which  is  used  in  calibration  of  sensor  values,  b,  The  histograms 
of  a  mixture  of  |  S )  and  |  T0)  states  used  to  calibrate  the  sensor  values.  If  we  choose  a  threshold 
Khresh  to  distinguish  between  |  S )  and  we  see  a  readout  fidelity  of  97%.  Purple:  fit  to  noisy 
distribution  including  Ti  decay  from  |  TQ)  to  |  S)  .  The  deduced  distribution  for  |  S)  (blue)  is  a 
Gaussian,  while  that  for  |  T0)  (red)  has  a  tail  due  to  Ti  decay. 
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3.3  Single  ST  Qubit 
Crosstalk  considerations 
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Fig.  3.3.4.  a,  A  schematic  of  the  energy  diagram  as  a  function  of  8  that  describes  the  two  qubits  and 
shows  the  regions  of  8  where  different  operations  are  carried  out.  b,  A  schematic  of  the  signal  from 
the  RF  charge  sensor  as  a  function  of  8  for  the  two  qubits.  This  signal  reflects  the  charge  distribution 
of  the  two  qubit  states.  For  large  positive  £  there  is  a  region  where  |  S )  and  |  T0)  have  the  same 
sensor  signal  (charge  distribution),  which  is  the  foundation  of  the  crosstalk-free  readout  scheme,  c,  A 
schematic  of  the  readout  scheme  that  eliminates  crosstalk.  First,  the  left  qubit  is  read  while  the  right 
qubit  is  "parked"  in  (0,2),  and  then  the  right  qubit  is  read  while  the  left  qubit  is  “parked”  in  (0,2).  d,  A 
two  dimensional  histogram  of  the  RF  sensor  responses  without  (left  panel)  and  with  (right  panel)  this 
crosstalk-free  readout  scheme,  e,  Histograms  of  sensor  values  without  (top)  and  with  (bottom)  the 
crosstalk-free  readout.  Without  the  crosstalk-free  readout  the  sensor  signal  of  one  qubit  depends  on 
the  state  of  the  other  qubit. 
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3.3  Single  ST  Qubit  -  Readout  Fidelity 


-0.015  -0.01  -0.005  0  0.005  0.01 


Probability  of 
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Fig.  3.3.5  We  prepare  |S>  many  times  and  measure  it.  We  also  prepare  |!>  many  times  (using 
prepared  hyperfine  gradient  dBz)  and  measure.  We  fit  the  sensor  signal  to  a  two  peaked  Gaussian, 
allowing  for  !1  (measured  independently)  decay,  and  allowing  for  imperfect  preparation  of  |!>  (i.e. 
some  of  the  states  that  we  initially  assumed  to  be  |!>  were  actually  |S>  for  instance  because  of  a  bad 
dBz  pulse).  We  then  pick  the  sensor  threshold  that  optimally  discriminates  between  a  hypothetical 
population  of  perfect  singlets  from  perfect  triplets  and  calculate  error  probabilities  from  the  measured 
distributions. 
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3.3  Single  ST  Qubit 

Single  ST  qubit  state  tomography 


a  is)  d  is) 


Fig.  3.3.6  Calibrated  versus  uncalibrated  state  tomography:  a,  Data  taken  to  calibrate  the  tomography 
shows  ripples  in  the  length  of  the  Bloch  vector  if  we  assume  that  the  tomography  projects  the  quantum 
state  on  to  Cartesian  axes  (inset),  b-c,  The  paths  around  the  Bloch  sphere  for  the  different  evolutions 
that  are  used  for  tomography  calibration.  If  the  tomography  is  assumed  to  project  on  to  the  Cartesian 
axes  there  are  points  that  lay  outside  the  Bloch  sphere,  and  the  pure  states  are  not  at  the  north  pole, 
which  is  indicative  of  flawed  state  tomography,  d,  The  ripples  in  the  length  of  the  Bloch  vector  are 
diminished  (compared  to  panel  a)  if  the  axes  deduced  from  state  tomography  (inset)  are  used,  e-f,  The 
paths  around  the  Bloch  sphere  for  the  different  evolutions  that  are  used  for  state  tomography.  When 
the  correct  axes  are  used,  all  the  points  lie  inside  the  Bloch  sphere  and  the  pure  are  at  the  north  pole. 
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3.3  Single  ST  Qubit 

Single  qubit  process  tomography 


Prepare  and  Measure  Initial  State 


-H 

•  •  • 

Prepare  Initial  Tomographic 

State  Measurement 


Run  Process  and  Measure  Result 


Fig.  3.3.7  Process  tomography  of  a  single  ST  qubit  uses  6  initial  states.  Processes  are 
repeated  many  times,  then  Chi  matrices  are  calculated  which  maximize  the  likelihood  of  an 
observed  process. 


R  tJ2 
2 


Fidelity  0.954 


Fig.  3.3.8  For  tuning  up  a  Clifford  set  we  use  a  set  of  axes  that  is  rotated  slightly  compared 
to  the  ST,  up-down,  axes  because  rotations  around  the  z  and  x  axes  are  challenging.  Shown  is 
the  7i/2  process  tomography,  tuned  around  the  new  z-axis,  which  is  the  axis  [.26,  .03,  .97] 
instead  of  [0,0,1].  Using  the  procedure  above,  we  calibrate  a  7i/2  rotation  around  this  axis  and 
find  a  process  fidelity  of  95%.  Shown  is  the  Chi  matrix  for  this  process  in  the  normal  S-T, 
up-do wn/down-up,  basis. 
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3.3  Single  ST  Qubit 
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Fig.  3.3.9  Process  Fidelity  for  Clifford  gate  set  indicating  the  axes  used  for  rotation.  J  is  7i/2 
rotation  around  the  designated  z  axis;  H  is  Hadamard  gate.  Note  example  fidelities  are  in  the 
range  >  0.94. 
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3.4  Two  ST  Operations 
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Fig.  3.4.1  The  two-qubit  pulse  sequence  used  to  entangle  the  qubits:  initialize  each  qubit  in  the  singlet 
state;  perform  a  tt/2  rotation  around  the  x  axis;  allow  the  qubits  to  evolve  under  exchange  for  a  time  t/2; 
perform  a  tt  rotation  around  the  x  axis,  thereby  decoupling  the  qubits  from  the  environment  but  not 
each  other;  evolve  under  exchange  for  t/2;  and  perform  state  tomography  to  determine  the  resulting 
density  matrix.  From  Ref.  [4]. 
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Fig.  3.4.2  (a)  The  elements  of  the  Pauli  set  of  the  measured  and  expected  states  for  t  =  140  ns 
(maximum  Bell  state  fidelity  0.72).  (b)  Maximum  fidelity  and  time  for  the  entangling  operation.  From 
Ref.  [4]. 
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Maximum  Fidelity  (%) 


3  S  Single  EO  Operations 


Fig.  3.5.1  Three  electrons  with  controlled  coupling  1-2  and  2-3  form  an  exchange  only  qubit  [6].  The 
two  control  axes  are  indicated  on  the  Bloch  sphere. 


Fig.  3.5.2  Readout  fidelity  of  EO  qubit.  Singlet  fidelity  is  0.96.  Triplet  fidelity  is  0.90,  lowered  because  of 
Ti  processes. 
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3  S  Single  EO  Operations 


state  tomography 


probabilities  in  non- 
orthogonal  basis 

Pi(p)  =  Tr[Eip\ 

P2(p)  =  Tr  [E2p] 
P3(p)=Tr[E3p] 
P4(p)  =  Tr[E4p ] 


five  density  matrices 

pi  =  ISiXSI 

P2  =  |5r)(5r| 

p3  =  Rr  Slurps i 

P4=  RtSr^RlSr 

P5  =  \mixed)(mixed\ 


Fidelities  for  the  Ef 

E\:  Fs  =  0.83,  FT  =  0.88,  V  =  0.71 

E2:  Fs  =  0.73,  FT  =  0.84,  V  =  0.57 

E3:  Fs  =  0.62,  FT  =  0.72,  V  =  0.34 

E4:  Fs  =  0.58,  FT  =  0.76,  V  =  0.34 

The  average  fidelity  of  the  triplet  outcomes  was  0.80. 

The  average  fidelity  of  the  singlet  outcomes  was  0.69. 

The  average  fidelity  of  all  measurements  was  0.75. 

The  average  visibility  was  0.49. 


Fig.  3.5.3  Pulse  sequences  for  the  EO  qubit  state  tomography. 
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